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Abstract 

High-temperature resummed perturbation theory is plagued by poor 
convergence properties. The problem appears for theories with bosonic field 
content such as QCD, QED or scalar theories. We calculate the pressure as 
well as other thermodynamic quantities at high temperature for a scalar one- 
component field theory, solving a three-loop 2PI effective action numerically 
without further approximations. We present a detailed comparison with 
the two-loop approximation. One observes a strongly improved convergence 
behavior as compared to perturbative approaches. The renormalization 
employed in this work extends previous prescriptions, and is sufficient to 
determine all counterterms required for the theory in the symmetric as well 
as the spontaneously broken phase. 
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1 Introduction and overview 



All information about the quantum theory can be obtained from the effective action, 
which is the generating functional for Green's functions. Typically, the (1PI) 
effective action T[<f>] is represented as a functional of the field expectation value or 
one-point function <ft only. In contrast, the so-called two-particle irreducible (2PI) 
effective action T[<p, D] is written as a functional of <fi and the connected two-point 
function D £Q. The latter provides an efficient description of quantum corrections 
in terms of resummed loop-diagrams. The different functional representations of 
the effective action are equivalent in the sense that they are generating functionals 
for Green's functions including all quantum/statistical fluctuations and agree by 
construction in the absence of sources. However, e.g. loop expansions of the 1PI 
effective action to a given order in the presence of the "background" field <fi differ in 
general from a loop expansion of T[<p, D] in the presence of (f) and D. 

This observation has been successfully used in nonequilibrium quantum field 
theory to resolve the problems of secularity and late-time universality El of 
perturbative approximations, which render the latter invalid even for arbitrarily 
small couplings |3]. Both the far-from-equilibrium behavior as well as the late-time 
thermal equilibrium results can be described from a loop expansion of the 2PI 
effective action 1 without further assumptions for scalar E] and fermionic 
theories j^j. Similar results have been obtained using a two-particle irreducible 
1/N expansion beyond leading order [HJ[71IE|. 

The same techniques can be applied directly in thermal equilibrium, where 
efficient formulations in Euclidean space-time become available. In contrast to the 
far-from-equilibrium case, there are various powerful approximation schemes known 
in thermal field theory. A prominent approach in equilibrium high-temperature 
field theory is the so-called "hard-thermal- loop" resummation jU]. However, explicit 
calculations of thermodynamic quantities such as pressure or entropy typically reveal 
a poor convergence except for extremely small couplings. An important example 
for this behavior concerns high-temperature gauge theories. Recent strong efforts 
to improve the convergence aim at connecting to available lattice QCD results, 
for which high temperatures are difficult to achieve. In order to find improved 
approximation schemes it is important to note that the problem is not specific to 
gauge field theories. Indeed it has been documented in the literature in great detail 
that problems of convergence of perturbative approaches at high temperature can 
already be studied in simple scalar theories. For recent reviews in this context see 
Ref. [TO]. 

A promising candidate for an improved convergence behavior is the loop or 
coupling expansion of the 2PI effective action. So far, thermodynamic quantities 
such as pressure or entropy have been mainly calculated to two-loop order. However, 
aspects of convergence can be sensefully discussed only beyond two-loop order since 
the one-loop high-temperature result corresponds to the free gas approximation. 
Efforts to calculate pressure beyond two loops include so-called approximately self- 
consistent approximations [H] 2 , as well as estimates based on further perturbative 

1 Loop approximations of the 2PI effective action are also called "<i>-derivable" . 

2 Cf. also Ref. J3| f° r a similar strategy in the context of Schwinger-Dyson equations. 
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Figure 1: Pressure as a function of the renormalized coupling, normalized to the ideal 
gas, i.e. one-loop pressure. Shown are the 2PI two- and three-loop results (left) for 
T = 2T with T = rriR(T ). The right figure shows the perturbative results to 
order g R and g\ along with the 2PI two-loop curve in the high-temperature limit 
for illustration of the problematic alternating behavior of perturbation theory (see 
text). 



expansions in the coupling and a variational mass parameter ^31 3 ■ These 
studies indicate already improved convergence properties. However, perturbatively 
motivated estimates as in Ref. suffer from the presence of nonrenormalizable, 
ultraviolet divergent contributions and the apparent breakdown of the approach 
beyond some value for the coupling. If one does not want to rely on these further 
assumptions, going beyond two-loop order requires the use of efficient numerical 
techniques. Such rigorous studies are important to get a decisive answer about the 
properties of 2PI expansions. As it turns out (cf. below) these problems appear as 
an artefact of the additional approximations employed and cannot be attributed to 
the 2PI loop expansion. 

In this work we calculate the pressure as well as other thermal quantities for a 
scalar g 2 (j) 4 field theory from a three-loop 2PI effective action numerically without 
further approximations. A detailed comparison with the two-loop approximation 
is presented. We observe a strongly improved convergence behavior as compared 
to perturbative approaches. This is exemplified in Fig. Q where the pressure is 
shown as a function of the renormalized coupling qr determined by the physical 
four- vertex. The left figure compares the two- and three-loop result normalized to 
the ideal gas pressure. For the employed high temperature T = 2m#(T ) the three- 
loop corrections to the pressure are rather small. Here mn(T ) is the temperature- 
dependent renormalized mass parameter or inverse correlation length and we have 
Tq = m^(To). For illustration we also show on the right of Fig. ^the perturbative 
results to order g\ and g\ along with the dominant 2PI two-loop result for the 
high-temperature limit. The problematic alternating behavior of the perturbative 



3 Cf. Ref. Jl] for a similar application to QED. 



series is not specific to the limit T ^> Tq and is characteristic for higher orders as 
well [13]. (Cf. Sees. O and EH ) 

We obtain the renormalized correlation functions or proper vertices as functions 
of temperature, building on a renormalization put forward in Refs. [TBI Ej (cf- also 
|18p. 4 In contrast to previous approaches, the renormalization employed in this 
work is formulated for the resummed 1PI effective action which is calculated from 
a given loop approximation of the 2PI effective action. The procedure is sufficient 
to determine all counterterms required for the theory in the symmetric as well as 
spontaneously broken phase. This extends previous prescriptions, which are not 
sufficiently general to renormalize the theory in the presence of a non- vanishing field. 
In particular, they do not renormalize all functional field-derivatives of the effective 
action even in the symmetric phase. The latter represent the proper vertices, which 
encode important information about the theory. This is discussed in Sec. 12. II taking 
into account three-loop corrections. The considered approximation represents a 
simple explicit example of a systematic renormalization scheme for 2PI effective 
actions, which will be given for general approximations in a separate publication j2JJ . 

To put our calculations in a more general context, we note that our results for 
the two- and three-loop approximations of the 2PI effective action are identical 
to the corresponding two- and three-loop approximations for n-particle irreducible 
(nPI) effective actions with arbitrary n > 3. The agreement up to the considered 
order is a consequence of an equivalence hierarchy for nPI effective actions [2*2"] . 
Their functional dependence takes into account the propagator as well as the proper 
three- vertex, four- vertex, . . . , n- vertex 01231 I2H 122] ■ Therefore, a loop expansion of 
the nPI effective action with n > 2 treats propagators and the respective vertices on 
the same footing, while a 2PI loop expansion singles out propagator resummation 
a priori. However, the difference to the 2PI results in the symmetric phase appears 
only at four- loop order [22], which is beyond the approximation employed in this 
work. As a consequence, the good convergence properties of the expansion, which 
we observe, may be attributed to the nPI loop expansion with arbitrary n > 2 rather 
than to the 2PI loop expansion. 5 

2 Renormalized thermodynamics 



We consider a quantum field theory with classical action 6 




4 For related studies see also Ref. ^5]. In contrast to the considered "non-local" 2PI 
resummation, these investigate renormalization for "local" resummations. The latter turn out 
to be problematic to describe the noncquilibrium late-time behavior of quantum fields 20 and will 
not be discussed here. 



5 An indication for this is that the quantitative description of the universal behavior near the 
second-order phase transition of this model goes beyond a 2PI loop expansion |25j . It requires 
taking into account vertex corrections that start with the 4PI effective action to four-loop order. 
The latter agrees with the most general nPI loop expansion to that order |22) . 

6 We employ a Minkowskian metric g^ v — diag(l, —1,-1,-1) in view of further possible 
applications, e.g. in out-of-equilibrium situations. 



where 4>(x) is a real scalar field with bare mass term m 2 and coupling g 2 . 
The normalization of the coupling is chosen for simple comparison with existing 
literature (cf. e.g. Ref. [TT]) in view of applications of these methods to QCD 
thermodynamics. We use the shorthand notation J = J Q dx° J d 3 x with 
temperature T. Following pP it is convenient to parametrize the temperature 
dependent 2PI effective action as 



V[<p, D] =S[<p] + ^Tr \11D- 1 + l -TiD Q \<p)D + T 2 [0, D] + const , (2.2) 

which expresses T in terms of the classical action S and correction terms including 
the function r 2 to which only two-particle irreducible diagrams contribute. Here the 
classical inverse propagator is given by iDQ 1 (x,y; </>) = 5 2 S[(j>]/ 8(j)(x)8(j)(y) . In the 
absence of external sources physical solutions require 



5T[(f>, D] 



5<f)(x 



5D(x,y) 



D=D(<f>) 



0, (2.3) 
0. (2.4) 



The 2PI effective action evaluated at D((p; x, y), i.e. for the ^-dependent solution of 
(I2.4JI . is identical to the 1PI effective action T[<f), D(<f))]. The effective action at the 
stationary point, T[(p Q , D((p )], corresponds to the logarithm of the partition function 
in the absence of sources . Therefore, in thermal equilibrium with temperature T 
(0o constant) the effective action is related to the pressure P by 

P=^r[0 O)J D(0o)], (2.5) 

where L 3 = J d 3 x denotes the spatial volume and the constant in Eq. ()2.2|) is chosen 
such that the pressure vanishes at zero temperature. Entropy density S and energy 
density £ are given by 

dP 

S = W , (2.6) 
£ = -P + TS = T 2 ^-(^] . (2.7) 



dT \T, 

We recall that all the physical information is contained in the effective action at the 
stationary point T[<po, D(0q)] and its changes with respect to variations in the field 
<p evaluated at <\> = (f) . For instance, the connected two-point function, T^ 2 \ and 
proper four-point function, r^ 4 ), are given by 

r»(,, y) , , (2.8) 

All information about the quantum theory can therefore be conveniently obtained 
from r[0, D(<f))} by functional differentiation. In particular, D{<fi)\ evaluated for 
constant field encodes the effective potential. 



2.1 Renormalization 



The effective action T[(f), D ((/>)] is defined in a standard way by suitable 
regularization, as e.g. lattice regularization or dimensional regularization, and 
renormalization conditions. We employ renormalization conditions for the two-point 
function ()2.8|) and four-point function (|2.9|) . which in Fourier space read: 

ZT^ 2 \p 2 )\ p=0 = -m 2 R , (2.10) 

Z ^ r(2)(p2)l -° = - 1 ' (2 ' n) 
Z 2 r^( Pl ,p 2 ,p 3 )\ Pl=P2=m=0 = -A\g 2 R , (2.12) 

with the wave function renormalization Z. Here the renormalized mass parameter 
m R corresponds to the inverse correlation length. The physical four-vertex at zero 
momentum is given by g R . Without loss of generality we use throughout this work 
renormalization conditions for <fio = 0. 

We emphasize that at non-zero temperature the mass parameter as well as the 
four- vertex are temperature dependent. The value of the mass and coupling at a 
given temperature and momentum scale uniquely determines the theory. This can 
then be used to calculate properties at some other temperature. We will typically 
define the theory by giving renormalization conditions at zero temperature as well 
as some non-zero temperature T with m fi = m R (T ) and g R = g R (T ). 



2.1.1 2PI renormalization scheme to order g R 

In order to impose the renormalization conditions (j2.1(Jj) - (j2.12j) one first has to 
calculate the solution for the two-point field D((j>) for = 0, which encodes the 
resummation and which is obtained from the stationarity condition for the 2PI 
effective action (|2.4|) . To achieve this we will for the most part follow the lines of 
the renormalization procedure described in Ref . . We comment on the additional 
implications, which arise from imposing (j2.10j) - (|2.12j) . below. For further details we 
refer to Ref. [21]. Here we will explicitly demonstrate from the numerical solution 
of the three-loop approximation that the renormalized quantities are insensitive to 
the change of the (lattice) regularization. 

The renormalized field is <p R = Z~ l / 2 <p. It is convenient to introduce the 
counterterms relating the bare and renormalized variables in a standard way with 

Zm 2 = m\ + 5m 2 , Z 2 g 2 = g\ + Sg 2 , SZ = Z - 1 , (2.13) 

and we write 

D{(j>) = ZD R {<j) R ) . (2.14) 
In terms of the renormalized quantities the classical action (J2.1|) reads 

S = f x {^d,4> R d^ R -\m 2 R <t> 2 R -gl^ R 

+ l -5Zd,ci )R d^ R - Um 2 <j>l - 5g 2 ^ R \ . (2.15) 



+ ® X ® 



Figure 2: Diagrammatic representation of the mass and field strength counterterms. 
The cross in the left graph denotes indistinctly the mass counterterm bm\ or field 
strength counterterm 8Z\ appearing in Eq. (j2.16|) . while the one in the right graph 
represents the 5m 2 and 5Z 2 insertions in Eq. (|2.15|) . The line of the closed loop 
represents D R , while a circled cross denotes (f) R . 



Similarly, one can write for the one-loop part TrlnD -1 = TrlnD^ up to an 
irrelevant, temperature independent constant. The next term of Eq. (|2.2j) can be 
written as 

7^r D (j 1 ((f))D((J)) = ~J (n x +mj t + 6Z 1 D x + 8ml)D R (x,y-,<p R )\ x=y 

-6 {g 2 R + 5gl) f <p 2 R (x)D R (x, x; <f) R ) . (2.16) 

J X 

Here SZx, 5m\ and 5gf denote the same counterterms as introduced in (|2.13jh 
however, approximated to the given order. To express T 2 in terms of renormalized 
quantities it is useful to note the identity 

r 2 [<p,D(<p)]\ g2 = T 2 [4> R ,D R {4> R )]\ gW , (2.17) 

which simply follows from the standard relation between the number of vertices, lines 
and fields by counting factors of Z. Therefore, one can replace in T 2 the bare field and 
propagator by the renormalized ones if one replaces bare by renormalized vertices as 
well. We emphasize that mass and wavefunction renormalization counterterms, SZ 
and 5m 2 , do not appear explicitly in T 2 . This can be understood from the fact that 
the only two-particle irreducible diagrams with mass and field strength insertions 
are those displayed in Fig. |21 The counterterms in the classical action (J2.15)) . in the 
one-loop term ()2.16|) and beyond one-loop contained in T 2 have to be calculated for 
a given approximation of T 2 . Here we consider the 2PI effective action to order g R 
with 



r 2 [4> R ,D R {4> R )} = - 




+zl2g R [ D%(x, y- <j> R ) - 35g 2 2 I D 2 R (x, x; fa) , (2.18) 

J xy J x 

where the last term contains the respective coupling counterterm at two-loop. There 
are no three-loop counterterms since the divergences arising from the three-loop 
contribution in ()2.18|) are taken into account by the lower counterterms fT7\ I2T] . 
The coupling counterterms are displayed diagrammatically in Fig. El 

One first has to calculate the solution D R (<ft R ) obtained from the stationarity 
condition ()2.4j) for the 2PI effective action. For this one has to impose the same 




+ 



Figure 3: Diagrammatic representation of the coupling counterterms. Here the 
vertices denote 5g%, 5gf and 5g 2 , respectively, which appear in Eqs. (j2.18|) . (j2.16j) 
and (l2~ToT) . 



renormalization condition as for the propagator (j2.10|) in Fourier space: 

iD R l (p 2 ; 0/?)| P =o,fe=o = -m% , (2.19) 
for given finite renormalized "four-point" field 7 

V R (x,y;z,w) = x ,,,,,, , • (2-20) 
b(f) R {z)(f) R (w) <j> R =o 

For the above approximation we note the identity 



S(j) R (x)(j) R (y) 

for 



= iDl L {x,y^ R )\^ R=Q (2.21) 

<t>R=0 



5Z = 5Z 1 , 5m 2 = 5ml > Sgl = 5gj, (2.22) 

such that (j2.19|) for D R is trivially fulfilled. In contrast to the exact theory, for 
the 2PI effective action to order g R a similar identity does not connect the proper 
four- vertex with V R . 8 Here the respective condition for the four-point field V R in 
Fourier space reads 

V R (p 1 ,p 2 ,P3)\ Pl = P 2= P 3=o = -&9r- (2-23) 

Note that this has to be the same as for the four- vertex (j2.12|) . For the universality 
class of the 4 theory there are only two independent input parameters, which we 
take to be m R and g R , and for the exact theory V R and the four-vertex agree. 
The renormalization conditions (j2.1(J)) - (J2.12|) for the propagator and four-vertex, 
together with the 2PI scheme ()2.19|) - ()2.23|) provides an efficient fixing of all the above 
counterterms. In particular, it can be very conveniently implemented numerically, 
which turns out to be crucial for calculations beyond order g 2 R . 

The apparent insensitivity 9 of the renormalized quantities with respect to 
changes in the (lattice) regularization is demonstrated for the temperature 

7 Below in Eqs. (|2.26(1 and Ij2.27|l we will see that for the present approximation this is the 
four-point field that solves the Bethe-Salpeter type equation discussed in Refs. [TCIIT7| . 

8 We emphasize that for more general approximations the equation (|2.21|) may only be valid 
up to higher order corrections as well. This is a typical property of self-consistent resummations, 
and it does not affect the renormalizability of the theory. In this case the proper renormalization 
procedure still involves, in particular, the conditions l|2.19|) and (|2.23|) . For a detailed discussion 
of these aspects see Ref. |2T] . 

9 Of course, we recover the "triviality" of 4 -theory such that the four-vertex vanishes in the 
continuum limit, which is discussed in Sects. 1^21 and EH 




Figure 4: The upper solid and dashed line show the results for Vr(T) as a function of 
the logarithm of the lattice cutoff for two temperatures T = 2T and T = 3T . Here 
Tq is fixed by T = m^(T ) and Vr(T ) = —1. For illustration we also show results 
for the same quantity (lower curves) if the renormalization is not done properly, by 
adding a "naive" three- loop coupling counterterm as explained in the text. (The 
"oscillation" is a lattice artefact arising from Fourier transformation of odd or even 
arrays.) 

dependent four-point field Vr{T) in Fig. |U The upper solid and dashed line show 
the results for Vr(T) as a function of the logarithm of the lattice cutoff for two 
temperatures T = 2T and T = 3T . Here T is fixed by T = m/j(T ) and the 
lattice volume is kept constant with Li71r(Tq) = 2. For illustration we also show 
the behavior of Vr(T) if the renormalization is not done properly. This is achieved 
by replacing in the three-loop contribution of (j2.18J) g R by the "bare" coupling 
(gji + Sg^) 2 and taking 5g1 = 5gf. The induced strong logarithmic cutoff dependence 
of Vr(T) can be observed from the respective results displayed in the lower curves 
of Fig. E| 

We emphasize that the current approximation (j2.18|) for the 2PI effective action 
can only be expected to be valid for sufficiently small (pR <C niR/gR. If the latter 
is not fulfilled there are additional 0(g R ) contributions at three- loop ~ g^R an d 
~ gR<pR- The approximation should therefore not be used to study the theory in 
the spontaneously broken phase or near the critical temperature of the second-order 
phase transition. For studies of the latter using the 2PI effective action see Ref . . 
Here we will investigate how the 2PI effective action cures convergence problems of 
high-temperature perturbation theory, for which the approximation ()2.18|) provides 
the correct generating functional up to order g R corrections, including all three-loop 
corrections in the high-temperature phase. 



2.1.2 Renormalized equations for the two- and four-point functions 



^From the 2PI effective action to order g R we find with 5g\ = bg\ from (J2.22J) for 
the two-point function: 



iDZ 1 



r {x,y;<PR) 



[(l + SZ l )n x 



m. 



5m\ 



+ \2{g 2 R + 5gl) (D R (x, x; <f> R ) + <j) R {x))] 8{x - y) 
+i288g%D R {x, y; ^r)M^)Mv) + iWg%!? R {x, y; 4> R ) . 



(2.24) 

According to (|2.21|) this expression coincides with the one for the propagator 
5 2 T[(f)R, D R (<f) R )]/5<f) R (x)8(l) R (y) at (p R = 0. It is straightforward to verify this using 



5T[(f> R ,D R ((f> R )] _ ST[<j> R ,D R ] 



Scj) R (x) 



(2.25) 



which is valid since the variation of D R ((j) R ) with <p R does not contribute due to 
the stationarity condition (|2.4|) . The four-point field ()2.20|) in this approximation is 
given by 



V R (x,y;z,w) 



-2A(g R + 5gl)S(x - y)S(x - z)S(x - w) 



-12{ 9 2 r + 5 9 2 / Dr ^ X ' M 



S(f) R (z)S(j) R (w) 



<>R=0 



5(x - y) 



+i 



(24^) 2 



8(x — w)8(y — z) + S(y — w)5(x — z) 



S 2 D R (x,y;4>R) 



5<)) R {z)5(j) R {w 
Inserting the chain rule formula 

5 2 iD R (x,y;(j) R 



<i>R=0 



D 2 R {x,y\(t)R = 0) 



5cj) R (z)5(j) R (w) 



<i>R=0 



D R (x, u; <p R )V R (u, v; z, w)D R (v, y; <j> R ) 



(2.26) 



^2.27) 



>>R=0 



into Eq. (j2.26|) one arrives at the Bethe-Salpeter type equation for V R discussed 
in Refs. [THI H7j. Eqs. (|2.26j) and (|2.24|) form a closed set of equations for the 
determination of the counterterms SZi, 8m\ and 5gf. Together with (|2.22j) one 
observes that 5g 2 would be undetermined from these equations alone as employed in 
Refs. |16| IT7] . This counterterm is determined by taking into account the equation 
for the physical four-vertex, which is obtained from the above 2PI effective action 
as 



5(f) R (x)5(j) R (y)5(f) R (z)8(j) R (w) 



4>n=o 



-2A(g 2 R + 8g 2 )S(x - y)S(x - z)5{x - w) 



+V R (x, y; z, w) + V R (x, z; y, w) + V R {x, w; y, z) 



5iDx(x,y\<l> R ) 
5D R (z,w) 



5iD R {x,z;(j) R ) 
5D R (y,w) 



+ 



5iD R 1 (x 1 w\(j) R ) 



5D R (y,z) 



<P R =0 



[2.28) 



where from (|2.24j) one uses the relation 



5iD R (x,y;(f> R ) 
SD R (z,w) 



<f>R=0 



12(g 2 R + 8gl)8{x - y)5(x - z)5(x - w) 



,(24<7i) 2 



D 



2 R (x, y- R = 0)6{x - z)5{y - w) .(2.29) 



2 



We emphasize that the counterterm 5g 2 plays a crucial role in the broken phase, since 
it is always multiplied by the field 0o an d hence it is essential for the determination 
of the effective potential. It is also required in the symmetric phase, in particular, 
when one calculates the thermal coupling using Eq. (|2.28p . 

2.2 Analytical example: 2PI two- loop order 

It is instructive to consider first the 2PI effective action to two-loop order for which 
much can be discussed analytically. 10 In this case Z — 1 for the scalar theory and 
the renormalized vacuum mass uir = m R (T = 0) of Eq. ()2.10|) or, equivalently, of 
Eq. ()2.19|) is given to this order by 



where m 2 = m R + 5m\ and g 2 = g R + bg\ is the zero temperature coupling 
Qr = 9r(T = 0) and we have used 1)2.22)1 . Here we have employed dimensional 
regularization and evaluated the integral in d = 4 — e for Euclidean momenta k. 
The bare coupling in the action (j2.1j) has been rescaled accordingly, g 2 — > fi € g 2 , and 
is dimensionless; Jx 2 = 47re~ 7B /i 2 and 7^ denotes Euler's constant. Below we will 
employ a lattice regularization for comparison and to go beyond two-loop order. 

Similarly, the zero-temperature four-point function resulting from Eq. ([2.2 3)) for 
the 2PI effective action to order g 2 R is given by 



with ()2.22)) . We emphasize that the same equation is obtained starting from the 
renormalization condition for the proper four- vertex ()2.12)) with 5g 2 = 35 g\. One 
observes that all counterterms are uniquely fixed by the renormalization procedure 
put forward in the previous sections. 

At non-zero temperature one obtains renormalized equations for the thermal 
mass and coupling in terms of the vacuum parameters: 




(2.30) 




(2.31) 



m 2 R (T) 



m 2 R - 24^ (T, m 



Cf. also |261 ITT] and references therein. 



+ 3 S ( m « (r) iog ^ - m ^ (T) + m *) ' (2 - 32) 



1 



Vr(T) A\g 



T -P>>(T,m 2 R (T)), (2.33) 



R 



9l{T) ~ l/^ + 4!^(T,m^(T)) (2 ' 34) 
where P denotes the pressure of the free gas: 

P (T, m 2 R (T)) = -T J J^L In (l - e^^) . (2.35) 

with cj p (T) = ^p 2 + m|(T) and P^{T,m R {T)) and P^'{T,m 2 R {T)) denote its first 
and second derivative with respect to m R (T). 
The pressure to two-loop order is given by 



P 2 = P (T,m R (T)) - - m 2 R (T)P>(T,m 2 R (T)) + 

-M4 + i) K(T)_m|)m| - (2 - 36) 

We have used these analytical formulae to check the (two-loop) numerics for the 
continuum and thermodynamic limit. We have also checked, that for the employed 
parameters the zero temperature limit is already reached for T < m R (T )/8. We 
use this below to numerically estimate the observables in vacuum. 

The scalar 4 -theory in 3 + 1 dimensions is non-interacting if it is considered 
as a fundamental theory valid for arbitrarily high momentum scales. This is the 
so-called "triviality" of 4 -theory [27]. However, if the theory is considered as 
a low-momentum effective theory with a physical highest momentum then the 
renormalized coupling gR can be non-zero. A non-vanishing zero-momentum four- 
vertex qr requires both a highest momentum cutoff and an "infrared cutoff" such 
as a non-zero mass term m R . This can be conveniently discussed by introducing a 
scale A with 

= \ + (2-37) 
In terms of this scale the renormalized coupling (four-vertex) reads 

g 2 = — — . 2.38 
yR 31n^- V ; 

As a consequence, the zero-temperature coupling g 2 R vanishes for ttlr/A — > 0, i.e. by 
either sending A — > oo for fixed tur or by taking tur — > for fixed scale A. For 
given A/rriR the renormalized coupling takes on its finite maximum value (|2.38|) for 
a diverging g due to the Landau pole in the equation for the coupling counterterm 
Sg 2 = g 2 — g 2 R obtained from (|2.31j) (cf. also the corresponding Fig. [7J for the 2PI 
three-loop solution). 

The high-temperature result can be obtained with non-zero tur for the limit 
T ^> rriR, which is conventionally dubbed the massless limit in the literature. To 



calculate m R (T) it is then useful to rewrite the equation (|2.32j) in terms of the 
renormalized four-point function at non-zero temperature, g R (T), for which the 
limit m R — > can be taken directly 11 (cf. also the discussions in Ref. ^U])- Here the 
two-loop pressure (|2.36|) in the limit m R — > becomes 



P 2 (T) = P (T) + \m 2 R (T)P^T, m|(T)) + ^± , (2.39) 



using 

m 2 3 A 

lim = lim — -ml In = . (2.40) 

m R ^0 gf, m R ^o 2ir 2 m R 

It is instructive to compare at this point to perturbation theory since its 
characteristic problems are apparent already at low orders. In Sec. 13.21 we will 
discuss aspects of the higher order behavior. To make link with other schemes we 
note that the MS renormalized "running" V R {ji) is given by 

4! 1 3 , m R 1 3 , A . 

--T + 7^1n-^ = -7 + 7^-, (2.41) 



V R {p) g\ 2vr2 fl g * 2n 2 e 
which can be written explicitly as 

- 7^00 = S • (2-42) 

4! l + &ln^ 

One observes that V R (p) corresponds to the zero- momentum four- vertex g R for the 
choice of fi = m R . 

In perturbation theory the massless or high-temperature limit has been 
particularly extensively studied in the literature. Up to order g R corrections the 
weak coupling expansion of the pressure for T 3> m R can be obtained starting from 
f|2.36|) or ()2.39|) using the high-temperature expansion of the gap equation ()2.32|) : 

ml(T) ~ g 2 R T 2 - —g R Tm R (T) , (2.43) 

7T 



where we neglect perturbative terms of order g R and higher. The first term comes 

(27T) 3 



from l2g R J 7^rn(p)/|p| = g R T 2 and the second term employs 



(2tt) 3 V ^ P |p| J ' J (2tt) 3 \p 2 + m 2 R (T) p 2 J 4vr 

with n(u p ) = (exp(oUp/T) — 1) _1 and where it was used that the momentum integral 
is dominated by momenta ~ g R T with n(u) p ) ~ T/u p . Using these approximations 

11 Note that at the critical temperature T c of the second order phase transition, determined by 
fnR{T c ) — 0, the temperature dependent coupling vanishes as well, i.e. g\(T) — > 0. However, 
the effective coupling of the dimensionally reduced theory g\{T c )T c /mn(T c ) remains non-zero and 
finite. We emphasize again that in order to quantitatively describe the universal behavior of the 
theory near T c requires to go beyond the present approximation (cf. e.g. the discussion using the 
l/N expansion of the 2PI effective action to next-to-leading order employed in Ref. and 
footnote 0) 



in the 2PI two-loop pressure (|2.39J1 . the high-temperature perturbative result is 
obtained as 

T 4 90 48 12tt ' v ' ' 

up to terms of order g R [10J. The entropy and energy density, respectively, are 
then given by 5 pert (T) ~ 4P pert (T)/T and £ pert (^) - 3P pcrt (T). The perturbative 
results to order g R and are displayed in Fig.[T]along with the 2PI two-loop result. 
The alternating behavior and poor convergence of the perturbative series observed 
from the low orders in (|2.44j) manifest itself also at higher orders of g R jT3]. This 
problem of perturbation theory is not specific to the limit T ^> m R but can also 
be observed for lower T/m R . For instance, for T <C m R the perturbative (positive) 
order g R contribution to the pressure is found to dominate the (negative) order g R 
contribution for not too small coupling [To] . 



3 Numerical calculation to three-loop order 

In order to obtain the results for the 2PI effective action to three-loop order without 
further approximations, we use numerical (lattice) techniques. The employed 
hypercubic lattice discretization provides a regularization scheme. It turns out to 
be convenient to carry out the calculation using the unrenormalized propagator D 
and four-point field V, where V = Z~ 2 Vr defined with Eq. ()2.20|) . On the lattice 
we consider Euclidean space-time and in Fourier space we denote the Euclidean 
two-point field by D k and four-point field by V k = V(k, — k, 0, 0), without loss of 
generality. Here D^ Q = Z~ l m R) dD^ 1 /dk 2 \ k=0 = Z~ x and t4=o — 24Z _2 g| ? are 
the Euclidean equivalents of the above renormalization conditions. Following the 
procedure of Sec. 12. 1\ for the renormalization one needs with ()2.22|) to know g 2 , 
Sg 2 , m 2 and Z which are functions of the lattice cutoff. The errors arising from 
subtracting cutoff dependent quantities such as the quadratic mass contribution 
of the setting sun diagram remain under control, since the numerical value of the 
unrenormalized quantities as well as the physical values simultaneously fit into a 
double precision variable for the employed cutoff values. 

On the lattice there is only a subgroup of the rotation symmetry generated by the 
permutations of p x ,P y ,p z and the reflections p x <-> — p x etc., which entails a reduction 
of independent lattice sites by a factor of 48. For three space dimensions, using 
periodic boundary conditions, a lattice of typical linear size 32a requires iV s = 969 
independent sites. In the fourth dimension we use a different lattice spacing (a t ) and 
lattice size (typically N t = 128) so the rotation symmetry cannot be extended. We 
implement two routines that make use of the symmetry features in order to have 
acceptable performance, the addup function and a four dimensional fast Fourier 
transformation defined as: 

addup(^) = ^V^, (3.1) 
fft(H) k = a 3 a t J2H x e ikx , (3.2) 

X 



invfft(i3% = —Y^H k e- lkx , (3.3) 
L4 k 

where L 4 = N 3 N t a 3 a t and H is a given lattice field in Fourier space. 12 Fields with 
a tilde (H) are in coordinate space. Here E& (and £ x ) denotes a summation over all 
lattice four-momenta (or coordinates) respecting the appropriate symmetry weights. 



3.1 Two-loop solution 

We first consider the two-loop lattice calculation, which can be compared with the 
analytic discussion of Sec. 12.21 In this case D^ 1 = k 2 + m 2 R and the respective 
coupling and mass equations on the lattice read 

l/g 2 = l/<£ -12addup(D 2 ), (3.4) 
m 2 = m \ - I2g 2 addup(D) . (3.5) 

We solve for g 2 and m 2 for given renormalized parameters as a function of the 
lattice cutoff. With these bare parameters we calculate the pressure, thermal mass 
and coupling at various temperatures, using the same lattice spacing. The results 
are given below along with the order g R results. The variation of the temperature 
is implemented by changing the lattice size according to T = l/(N t a t ). 

The pressure is calculated by evaluating the 2PI effective action to this order 
using Eq. ([2.5)1 : 



P2 = addup 



^log£> + 3g 2 D&ddupD 



(3.6) 



This expression suffers from a temperature independent quartic divergence. To 
subtract this, we always measure pressure differences, by calculating the pressure 
at zero temperature as well. We carry out the subtraction before calling the spatial 
part of the overall addup function. The sum over the fourth dimension has to be 
done beforehand, since the respective lattice sizes are not the same for different 
temperatures. 



3.2 Three-loop solution 

We obtain the renormalized physical results in two steps: At zero temperature or 
a given temperature T we numerically solve the equation ()2.24|) for the two-point 
field and the equation ()2.26|) for the four-point field with the conditions (j2.19|) and 
(12.23)1 simultaneously. This way we obtain the counterterms Sgf, 5Zi and 5mi, using 
()2.22j) . The counterterm for the coupling in the classical action, 5g 2 , can then be 
obtained from Eq. (|2.28|) . At a different temperature T we use the counterterms 
obtained in the first step to evaluate physical quantities. We solve first Eq. ()2.24)) 
and get the thermal propagator, the thermal mass and the pressure. Then by solving 
the corresponding equation ([2.26)1 we obtain the four-point function Vfc(T). Finally, 
the thermal coupling is obtained from Eq. ()2.28[) using the previously obtained value 
of«V. 

12 For x — the values of Eqs. i|3.1fl and i|3.3fl are identical. Both routines are implemented 
separately to reduce computational costs. 



In the following we describe the numerical implementation of the simultaneous 
iterative solution of the propagator and coupling equation in more detail. The 
renormalized mass, uir, is extracted from the low-momentum form of the inverse 
propagator D^ 1 = Z^ 1 (k 2 + m 2 R + C(/c 4 )). The renormalized coupling, g 2 R is 
obtained from the respective vertex equation for V k to this order: 

ya = ^ 2 -^E^- 24 ^ 4 ^ + ^^Ew^> (3-7) 

4 q 4 q 

with F k = L^YjqDqDk-q- We note that both the integrand and the kernel F k from 
the sunset diagram can be calculated by convolution, which is simple if the f f t and 
invf f t routines are provided. The numerical implementation of the vertex equation 
can be summarized as follows: 

F x = (invfft(D))* 

F k = fft(F) fc 

begin loop (iterations) 

H x = invift(D 2 V) x 

I k = fft(FH) k 

= g 2 - g 2 H x=0 /2 - 2Ag R Z' A F k + I2g% Z~H k 

end loop 

Similarly, the propagator equation can be implemented by convolutions: 
begin loop (iterations) 

D = invfft(D) 

= D -l + 12^2 addup p) _ 9 Qg R Z^fft(D 3 ) k 

end loop 

with DqI = k 2 + m 2 . We emphasize that the simple iterations do not converge 
well. The iterated propagator and four-point function oscillates around a physically 
senseful value and sometimes this oscillation is damped slower than the round- 
off errors accumulate. The alternating behavior at each order in the iteration 
is exemplified for the zero-momentum propagator in Fig. A more efficient 
iterative solution procedure involves a convergence factor, which avoids/damps the 
alternating behavior: 

For a choice of a = 0.1... 0.5 convergence was typically achieved after 20 . . . 50 
iterations by matching the exit criterium. This criterium is based on the sufficient 
smallness of the absolute change in D or V at all momenta. 

We note that the naive iterative solution of the equation for the propagator 
(with a = 1) reflects some aspects of a perturbative calculation. Starting from the 
classical propagator, each iteration adds a higher order contribution to D. At low 
orders the origin of the problematic oscillating behavior of a perturbative series in 
gu can be nicely observed from Fig. |S| Though higher orders in the iteration do not 
take into account all respective perturbative contributions, it is interesting that a 
convergence is obtained only after of order 100 (!) iterations. 

When fixing the renormalized theory we implement an outer loop for repeating 
the subsequent solutions of the vertex and the propagator equation. The bare 
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Figure 5: Shown is the behavior of the naive iterative solution of the equation for the 
propagator. Here each number corresponds to the value of D^ =0 after the respective 
number of iterations, starting from the classical propagator Z?o,&- The temperature 
is T = Amji in units of the zero-temperature mass rriR. One observes convergence 
only after of order 100 iterations. Iterative solution procedures can nevertheless be 
very efficiently applied with slight modifications (with a convergence factor of 0.5 
convergence is achieved after 14 steps, cf. the discussion in the text). 



coupling is tuned within the vertex equation iterations and the bare mass is adjusted 
after each iteration in the propagator equation. Then, using Eq. (|2.28j) . simple 
algebra gives 5g 2 . Calculating the thermodynamics at a given temperature using 
the previously obtained bare parameters is much simpler: first we solve the equation 
for the two-point field without any tuning of the parameters, then we obtain the 
thermal coupling from the vertex equation and no outer loop is required here. In 
addition to the two- loop contribution to the pressure, P 2 , there is an additional 
contribution from the basketball diagram at three-loop order such that the pressure 
is given by 



addup 



- log£> + 3g 2 Da.ddvpD - 36g%Z~ i Df f tD 3 



(3.9) 



with D = invf f tD. The subtraction of the divergences proceeds exactly in the 
same way as in the two-loop case: we are subtracting spatial lattice fields and we 
can only then carry out the spatial integral. 

Though this analysis is carried out using a lattice discretization, there is no 
principle obstacle exploiting the full rotation invariance of the continuum equations. 
This will be important in order to discuss the high-temperature behavior in 
applications of these techniques to more complex theories such as QCD. 
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Figure 6: The left figure shows the results for pressure at a lower temperature 
(T = ttir(O)) than Fig. El as a function of the renormalized zero-temperature 
coupling. The right figure shows the renormalized thermal coupling as a function of 
temperature. As for Fig. ^results are given employing two different renormalization 
scales. The shaded band indicates the renormalization scale dependence appearing 
at three-loop order. 



4 Discussion 



On the left of Fig. ^ the two- and three-loop results for the pressure are shown as 
a function of the renormalized coupling g R . For the employed high temperature 
T = 2T with T = m#(T ) we observe that the three-loop correction to the pressure 
is rather small. We note that here the three-loop thermal mass is larger than the 
two-loop mass, which drives the three-loop pressure below the two-loop value. 

In Fig. |Ul we show results at a lower temperature, which is taken to be equal 
to the vacuum mass, i.e. T = m R . We plot the pressure as a function of the 
vacuum coupling. This set of parameters has the feature that the thermal masses 
are within 1% equal in the two approximations for any coupling. In this case the 
three-loop curve is above the two-loop curve. Note that the two-loop contribution 
from the 2PI effective action is negative, while the three-loop correction is positive. 
This can be directly observed from the 2PI effective action before evaluation at 
the stationary point (|2.4|) for which iT^ oop [(f) = 0, D]T/L 3 = — 3g 2 (addupD) 2 and 
ir^ loop [0 = 0,D]T/L 3 = 12(^/Z) 4 addup(i)fftZ) 3 ). As a consequence, if the 
propagators do not change much from two-loop to three-loop order then the pressure 
will always be increased by the three-loop contribution. On the right of Fig. El we 
also show the renormalized coupling as a function of temperature normalized to rriR. 

For the exact theory the choice of a temperature scale for renormalization is 
irrelevant. However, for a given approximation the possible renormalization scale 
dependence of the results can be used as a check. The two-loop result is manifestly 
renormalization scale independent (cf. also Sec. 12. 2|) . The three-loop results of 
Figs. ^ and |H1 have been calculated first from renormalization at zero temperature 
and second from renormalization at a non-zero temperature T = ttir. The high- 
temperature renormalization conditions are chosen such that the results agree at 
zero temperature. At three loop this check is nontrivial because of the coupling 
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Figure 7: Physical coupling ranges for different lattice spacings employing the three- 
loop 2PI effective action. Here we show the coupling counterterms 5gf = bg\. The 
counterterms diverge at a certain value of the physical coupling. 

dependence in Eq. (I2.18|) with g R = Qr{T = 0) for the first calculation and 
9r — 9r{T = trr) for the second one. The variation of the thermal results of 
the two models give an idea about the scale dependence, which is indicated by 
the shaded bands in Fig. |Ul A similar analysis employing different renormalization 
scales has been also performed for the 2PI three-loop results displayed on the left of 
Fig.^ However, the difference between the results was hardly distinguishable in this 
high-temperature case. In contrast, the severe scale dependence of the perturbative 
calculations is well documented in the literature jlOj . 

Compared to the perturbative expansions for pressure (cf. Fig. ^ right), our 
results indicate a substantially improved convergence for the 2PI expansions. This 
holds even for rather strong couplings. In our lattice calculations we approach the 
Landau pole with our momentum cutoff, i.e. we explore almost the full range of 
couplings principally available. The range of renormalized couplings for various 
lattice regularizations is shown in Fig. For any momentum cutoff there is a 
highest value of g R for which there exists a bare coupling g 2 = g R + Sg 2 . (Cf. also 
the discussion for the two-loop effective action in Sec. I2.2I ) Concerning the lattice 
discretization one gains precision by reducing the lattice spacing, however, the range 
of physical couplings that can be realized shrinks. From the comparison to the 
two-loop analytic formulae, and from three-loop calculations with a series of small 
lattice spacings, we infer that one should use am^(T ) < 0.125 (for T < 2m^(T ), 
see Fig. HJ) . How this limits the available couplings to three- loop order is shown in 
Fig.0 

Very similar techniques as those employed here can be straightforwardly used 
for more complicated systems, as for instance fermionic or Yukawa theories (cf. also 
Ref. [B]). A more ambitious generalization is the application of these techniques to 
gauge field theories. While linear symmetries as realized in QED can be treated 
along similar lines, the generalization to nonabelian gauge theories is technically 
more involved and needs to be further investigated |21j . 
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